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ABSTRACT 

We report on a Suzaku observation of the newly discovered X-ray binary MAXI J1836-194. The source is 
found to be in the hard/intermediate spectral state and displays a clear and strong relativistically broadened 
iron emission line. We fit the spectra with a variety of phenomenological, as well as physically motivated disk 
reflection models, and find that the breadth and strength of the iron line is always characteristic of emission 
within a few gravitational radii around a black hole. This result is independent of the continuum used and 
strongly points toward the central object in MAXI J1836-194 being a stellar mass black hole rotating with a 
spin of a = 0.88 ± 0.03 (90% confidence). We discuss this result in the context of spectral state definitions, 
physical changes (or lack thereof) in the accretion disk and on the potential importance of the accretion disk 
corona in state transitions. 
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1. INTRODUCTION 

The X-ray spectra of X-ray binaries provide important clues 
on the nature of the compact objects and on the broad prop- 
erties of the accretion flow. In particular, the various re- 
flection features endemic to stellar mass black hole binaries 
in all active states have been successfully used to constrain 
the dimensionless spin parameter of various black holes in 
binary systems, as well as providing invaluable insight into 
the manner in which accreti on flow varies with mass accre- 
tion rate (e.g. lMilleiil2007h . By far the most prominent - 
and probably the most important - of these reflection features 
is the relativistic iron line appearing at approximately 6.4 to 
6.97 kcV depending to the ionisation state of the emitting ma- 
terial (e.g. Tanaka et al. 1995). Suzaku combines exceptional 
energy resolution below ss 10 kcV with broadband observa- 
tion, and, as such, is unique amongst other current X-ray satel- 
lites in the study of black hole transients. 

The observed spectrum often exhibit the presence of ther- 
mal emission, originating in an optically-thick accretion disk 
together with a hard component often referred to as the 
corona, and reflection features. There is a complex link be- 
twee n the disk, corona l hard X-rays and reflection emissions 
(e.g. iDone et alj|2007[). which manifests in various spectra l 
states (see lRemillard & McClintockll2006l and lBelloniir20I0h . 
Characterising the driving force between these state transi- 
tions is a fundamental challenge for both theoretical and ob- 
servational studies of accretion-flow properties. 

The prevailing paradigm requires that in quiescence (very 
low m), the inner accretion disk is fully replaced by an Ad- 
vectio n Dominated Accretion Flow (ADAF; e.g. lEsin et aU 
119971) . This has led to the idea that the transition between 
active states is a manifestation of changes in the innermost 
extent of the accretion disk, so that a transition from the disk- 
dominated high/soft state to a powerlaw-dominated low/hard 
state marks the point of the disk recession. The constant pres- 
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ence of radio jets in the low/hard state can also be associ- 
ated, albeit in a qualitative manner, with the truncation radius. 
Thus, it is clear that knowledge of the inner extent of the ac- 
cretion disk can have fundamental consequences to our under- 
standing of the nature of the accretion flow at low m as well as 
the connection between accretion disk, corona and radio jets. 
This radius can be determined via the study of both the con- 
tinuum emission from the accretion disk or by the reflected 
iron emission line, again making Suzaku, with its broad band 
coverage and high spectral resolution, ideal for this science. 
In fact, the advent of XMM-Newton and Suzaku has strongly 
challenged the paradigm that the accretion disk is truncated in 
the brig ht phases of the low/hard state in black hole binarie s 
(but see lDone & Diaz Trigol20Tot also see lMiller et al.l2010h . 

This challenge is exemplified by the recent Suzaku 
and XMM-Newt on observations of XTEJ1752-223 
dReis et al.ll2011bl) and flie 42 ks XMM-N ewton observation 
of XTE J1652^53 (lHiemsttaeLalJ|20li]). XTEJ1752-223 
was caught during the decay of its 2009 outburst in both 
the intermediate {Suzaku) and low/hard (XMM-Newton) 
spectral states. Interestingly, in both observations we 
found the presence of a strong, relativistic iron emission 
line which independently yielded strong constraints on 



the inner radius: 
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(90% confidence; iReis et alJ 1201 Ibl) for the intermediate 
and low/hard state respectively, as well as thermal disk 
components clearly following the L oc T'^ relation expected 
for geometrically-thin accretion disk, thus strongly ruling 
out disk truncation in either state. Similar results have been 
found for XTE J1652-453, where iHiemstra et alJ JMTh find 
the disk to be at w 4rg in the hard/intermediate state, and 
for GX 339-4 where th e disk doe s not appear to truncate 
until at least 10~^L rw^ (iMilleret al . 2006; Reis et al. 200i 
iTomsick et a l."2009'; Wilk inson & Uttlevll2009l) . 

In order to determine whether these sources are anomalous, 
or if state transitions are really not linked with a recession of 
the innermost extent of an optically-thick, geometrically-thin 
accretion disk, we urgently need observations at low fractions 
of the Eddington limit, down to 10^** and below. However, it 
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is clear tha t at in some phases of the low/hard state, at least as 
defined by iMcClintock et aTl (l2006b . the disk does not appear 
to truncate beyond the radius of the Innermost Stable Circu- 
lar Orbit (ISCO). Whether this is due to a different phase of 
the low/hard state - i.e. an "ADAF state" - or whether the 
disk only truly begin to truncate at much lower rii remains 
to be seen. Either way, it is clear that the disk plays an im- 
portant role in state transitions and even jet creation, be it by 
physically truncating and allowing for the existence of an in- 
ner ADAF zone or by dissipating less gravitational energy and 
allowing for a more powerful accretion disk corona. 

The shape of the iron line is determined by the relative 
depth of the disk within the potential well of the black hole 
and, as such, conveys information on its spin. Understanding 
the role of black hole spin (a = cJ/GAP, -1 < a < 1) 
in shaping accretion flows onto and jets from black holes is 
an important goal, having strong repercussion in all areas of 
astronomy. Stellar-mass black holes, for example, are likely 
to gain most of their angular momentum during birth, and 
their spin is a consequence of the supernovae t hat results in 
the cr eation of the central black hole (see e.g. iMiller et"aI1 
1201 II) . Knowing the spin distribution for these objects thus 
provides a window into the nature of one of the most power- 
ful explosions in the Universe. At present, we have approx- 
imately a dozen spin me asurements m ade by the use of the 
relativistic iron fines (e.g. iMilleret af1 l2008. 2009; Rei s et al.l 
120081 12009al 1201 Ibh . and a handful obtained from the ther- 
mal disk continuum dMcClintock et a} . 2006; Gou et al. 20091 
iSteiner. Reis. McClintock et al.H2011l) . However, in order to 
make any claim on the possible rol e of spin on, for example, 
radio jet power dPender et al.ll20ldl) . we need to increase our 
spin demographics. 

In this paper, we draw on the recent Suzaku TOO ob- 
servations of the nearly discovered black hole candidate 
MAXI J1836-194 to learn about the nature of the innermost 
accretion flow in this source, and to increase our black hole 
spin demographics. The following section summarises all the 
observations of the source and details the current Suzaku ob- 
servation. Section 3 begins by exploring some of the more 
phenomenological models used to explain the spectra of X- 
ray binaries and confirms the black hole nature of the cen- 
tral source. We conclude this section by using a fully self- 
consistent and physically motivated model to estimate the spin 
pai-ameter of the black hole in MAXI J1836-194. Section 4 
summaries our results and discusses the implications for cur- 
rent ideas of black hole state transitions and interpretations. 

2. OBSERVATION AND DATA REDUCTION 

MAXI J1836-194 was discov ered by the MAXI /GSC 
observatory on 2011 August 30 jNegoro et al.ll2011b . Its 
current evolution, at various wavelengths, have been reported 
in various ATels dK ennea et al. 2011; Cenko et al. 20111 
Sfi-ohrnaver & Smith 2011; Rau et al. 2011; Nakahira et at] 
201 Ibt iDel Santo et al., ,2011; .Trushkin et al., ,201 It 



Russell et al.l 1201 fl) with the latest report bv iRussell et al.l 
(120111) strondv suggesting that MAXIJ1836-194 is indeed 
a black hole X-ray binary based o n VLT mid-IR detections . 
Similar conclusions were made by IMiller- Jones et al.| ( 1201 II) 
based on EVLA radio detections. The 6 ks RXTE/PCA obser- 
vation of MAXI J 1 8 36- 1 94 reported by iSfi-ohmaver & Smithl 
( 1201 ll) had the spectrum described by an absorbed pow- 
erlaw with photon index of 1.84. The authors reported 
the presence of an iron emission line at 6.3 ± 0.2 kcV 
together with a smeared edge at 7.2 ± 0.2 kcV and alluded 
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Fig. 1.— SwiftTBAT lightcurve for MAXI J1836-194 in the 15- 
50koV energy range. The dashed horizontal lines shows the 
average rate and +/- 1 sigma standard deviation for this source 
since its discovery. The aiTow shows the time of the Suzaku 

observation reported in this paper Lightcurve available from: 

|http://heas arc.gsfc.nasa.gov/docs/swift/results/transients/weak/MAXIJ 1836-1 94/ 1 

to a reflection interpretation of these features. The quoted 
3-20 kcV flux of 9.8 x 10~^° erg cm~^ s^^ is remark- 
ably similar to that of XTE J 1752-223 in a si milar state 
(» 11 X IQ-i" erg cm-2 g-i; T = 1.83 ± 0.02; iReis et alj 
I2011H) . 

Figure [U shows the hard X-ray evolution of 
MAXIJ1836-194 as observed by the Swift/BAT Hard 
X-ray Transient Monitor provided by the Swift/BAT team. 
Suzaku observed the source on 201 1 September 14 for a total 
of 10.7 ks during a period where the 15-50 kcV flux showed 
a slight decrease. The three operating detec tors constitut- 
ing t he X-ray Imaging Spectrometer (XIS; iKovama et alj 
l2007h on-board of Suzaku were operated in the 1/4 win- 
dow "burst" mode with both front and back illuminated 
detectors in the 3x3 and 5x5 editing modes. Using the 
latest HEASOFT v6.11.1 software package we processed the 
unfiltered event files for each CCD following the Suzaku 
Data Reduction GuideQ. Due to the observation having been 
preformed in "burst" mode, we started by producing detailed 
good time intervals (GTIs) using the FTOOL XISTIME and 
setting the option "bstgti=yes". New attitude files were then 
created using the AEATTCOR scripfl ( lUchivama et al.ll2008h 
in order to correct for shift in the mean position of the source 
caused by the wobbling of the optical axis. The FTOOL 
XISCOORD was used to create new event files which were 
then further corrected by re-running the Suzaku pipeline with 
the latest calibration, as well as the associated screening 
criteria files. The good time intervals provided by the XIS 
team were also employed in all cases to exclude any possible 
telemetry saturations. XSELECT was used to extract spectral 
products from these event files. 

In order to estimate the level of pile up suffered by the data 
we used the script PILE_EST0 (lDavisll200lh to create a pileup 
map out of a Suzaku event data file. After experimenting 
with various extraction regions we choose to employ a box 
annulus region with a width of 240 pixels (~ 250") and a 
height of 290 pixels (~ 302") and an inner radius of 70 pix- 
els (^ 73"). This resulted in a maximum pileup fraction of 

^ http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/ 

http://space.mit.edu/cxc/software/suzaku/aeatt.html 
' http://space.mit.edu/cxc/software/suzaku/pest.htm 
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2 per cent. Background spectra were extracted from a cir- 
cular region having a radius of 100" elsewhere on the same 
chip. Individual ancillary response files (arfs) and redistri- 
bution matrix files (rmfs) were produced with the script XIS- 
RESlQ - which calls the tools XISRMFGEN and XISARFGEN 
- with the "medium" input. 

Finally, we combined the spectra and response files from 
the two front-illuminated instruments (XISO and XIS3) using 
the FTOOL ADDASCASPEC to increase signal-to-noise. The 
FTOOL GRPPHA was used to give at least 100 counts per 
spectral bin in a total of 512 energy channels. The nominal 
energy range covered by the XIS detectors is from ~ 0.2 — 
12 keV. However there are still calibration issues below ~ 
1 keV, which are especially severe in the burst clocking mode. 
Therefore, we do not consider data below 1.2 kcV or above 
lO keV, following the analyses of XTE J 1752-223 presented 
by iNakahira et alJ (1201 1 ah . The energy band of 1.6-2.4 kcV 
is also excluded to avoid large systematics uncertainties in the 
effective area near the silicon K and gold M edge. The smaller 
effective area, together with the fact that the out-of-time event 
rate is more significant in the BI instrument, as compared to 
the FI, means that at the energies considered here (i.e. above 
1 .2 keV) the BI data contain even larger uncertainties, and for 
this reason we did not use XIS 1 data in this paper 

We processed the The Hard X-ray Detector (HXD; 
iTakahashi et al.ll2b07l) with the standard criteria. The appro- 
priate response file (ae_hxd_pinxinomel 1_201 10601. rsp) for 
XlS-nominal pointing was downloadecQ and the data were 
reprocessed in accordance with the Suzaku Data Reduction 
Guide. As the non-X-ray background (NXB) file was yet to be 
created by the HXD team at the time of writing, we estimate 
the NXB by extracting the earth occulted data" (ELV< -Sfl 
Dead time corrections were applied with HXDDTCOR. The 
contribution from the Cosmic X-ray Backg round (CXB) was 
simulated using the form of iBoldti (119871) with appropriate 
normalisation for the XIS nominal pointing, resulting in a 
CXB rate of 0.019 comit s'^. The earth occulted NXB and 
CXB spectra were then combined using MATHPHA to give 
a total background spectrum, to which a 2 per cent system- 
atic uncertainty was added. The source spectrum was finally 
grouped to at least 100 counts per spectral bin. The PIN spec- 
trum is restricted to the 15.0-42.0 kcV energy range and fit 
simultaneously with the XIS data by adding a normalisation 
factor which is set to 1.16 in all fits with respect to that of 
the FI spectrum as recommenced by the Suzaku data analysis 
guide. To test the robustness of our result, in the following 
sections we also investigate the effect of allowing this cross- 
normalization to vary. All errors reported in this work are 
90 per cent confidence errors obtained by allowing all param- 
eters to vary, unless otherwise noted. 

3. DATA ANALYSES AND RESULTS 

EXPLORING PHENOMENOLOGICAL MODELS 

In order to compare the spectral properties of 
MAXIJ1836-194 with past work on other X-ray bina- 
ries, we begin by fitting the data with a simple combination of 

' http://suzaku.gsfc.nasa.gov/docs/suzaku/analysis/xisresp 
' http://www.astro.isas.ac.jp/suzaku/analysis/hxd/ 

^ This has now become available and we have checked the consistency in 
the results. The "tuned" background has a 15-45 kcV flux that is approxi- 
mately 3% less than the earth-occulted background, and is fully consistency 
within errors. The spectral shape is also fully consistent with one another and 
does not alter the results presented here. 




Energy (keV) 

Fig. 2.— Response unfolded i/Fu spectrum of MAXIJ 1836-1 94. The 
total, disk and powerlaw components are shown as cyan, blue and green solid 
lines respectively. Bottom: Data model ratio to a absorbed DISKBB plus 
POWERLAW fit ignoring the 4-7 kcV energy range. 



an ab sorbed powerlaw together with a DISKBB ( iMitsuda et al.l 
[1981 model. Fi gure 2 shows this fit with the 4-7 keV range 
ignored in order to best model the continuum. The total 
0.01-100 keVunabsorbed flux is - 7.8 x 10"^ erg cm^^ g-i 
of which approximately 42 percent is associated with the 
accretion disk. For com par ison with the standard work of 
iMcCUntock et al.l (120061) and lBeUonil (120101) . we quote the to- 
tal 2-20 kcV unabsorbed flux as ~ 1.5 x 10~^crg cm^^ s^^ 
and a disk fraction of w 0.26. Combined with a spectral index 
of 2.0, this observation of MAXI J1836-194 is consistent 
with having caught the source in the hard/intermediate 
spectral state. A further hint of this is also seen in Fig. 1, 
where it is clear that the source was observed during a time 
of slight decrease in the hard X-ray flux and suggests a short 
excursion away from the low/hard state. 

In Figure 2, the 4-7 keV energy range is shown above the 
continuum in order to highlight the presence of various fea- 
tures in this range. A possible explanation used by a number 
of authors to account for the residuals seen in Fig. 2 assumes 
that these features are a combination of a narrow iron Ka 
emission line and its associated absorption edge, which, when 
arising from the region around a black hole, suffers from a 
high degree of smeari ng and thus is better described by the 
SMEDGE model (e.g. lEbisawa et"anil994l) . A fit with such 
a smeared edge having an energy of 7.11 kcV (as expected 
from neutral iron) and a width of lOkcV was not able to ac- 
count for the residuals {x^/v = 630.9/328 Model la; Table 
1 and Figure 3), with broad residuals remaining both above 
and below the edge energy. Adding a narrow (ct = 1 eV) 
GAUSSIAN line at 6.4 keV did not improve the residuals in 
any way. A much better fit is indeed achieved when the ener- 
gies of both the GAUSSIAN line and that of the smeared edge 
are allowed to be free and the emission line is allowed to be 
broad {x^/v = 380.3/324; Model lb; Table 1). The neutral 
hydrogen density found in these models are also mildly con- 
sistent with the valu e of 2.0 ± 0.4 x 10^^ cm~^ presented by 
iKennea et alj ( 1201 lb based on a Swift/XRT observation. We 
will discuss possible reasons for the variation in A'h, as ob- 
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TABLE 1 

Phenomenological fits used to determine the robustness of the emission line. 





Model la 


Model lb 


Model 2 


Model 3a 


Model 3b 


Model 3c 


Model 3d 


Nil ( xl022 cm-2) 


0.08 ±0.02 


0.11 ±0.02 


05+"-"^ 
^■'■'°-0.02 


0.05 ±0.03 


0.05 ± 0.03 




< 0.05 


Eamedgn ( keV) 


7.11(f) 




— 


— 


— 


— 


— 




2.7 ±0.2 


3.25 ±0.3 


— 


— 


— 


— 


— 


r 


2.02 ±0.01 


2.05 ±0.01 


i.oau_g 


2.13 ±0.06 


2.17 ±0.03 


2.18 ± 0.02 


2.09 ±0.03 


fcTdi.fc (keV) 


0.445 ± 0.005 


0.433 ± 0.005 


0.474 ± 0.003 




0.30 ±0.01 


0.30 ±0.01 


0.30 ±0.01 


^^electron ( ke V) 


— 


— 


— 




10+2 


10+2 


10+s 




— 


— 


— 


0.358_q;o3 


0.39 ± 0.02 




n qo+O-OS 


^hard 


0.36 ±0.01 


0.38 ±0.01 


0.219 ± 0.001 


r. r,Y + 0.06 

—0.05 


Q 33+0.04 
— 0.03 


0.33 ± 0.03 


0.29 ±0.03 


A^disfe66(Xl03) 


4.4 ±0.3 


5.1 ±0.4 


3 54+0 03 












A^CompPs(Xl03) 


— 


— 




30t^ 


gi+16 
"^ — 14 




53tf 


EgaiLss (keV) 


— 


6.4-6.97 


— 




5.5 ± 0.1 




— 


Tjauss (eV) 










sool^r 






£^Iaor(keV) 






6.97-0.04 






f. c:o+0.14 


6.48+0.10 


Af,„e(xl0-3) 




2.6±0.7 


1,5 4+0.4 




2.8tO:« 


o 9+0.3 
■'•^-0.4 




W'iine(eV) 




300 ± 90 


360 ± 10 






270^20 


180t«« 








6.8 ±0.2 






n Q+0.3 


3.2 ±0.2 


8 (degrees) 






44 ± 1 






25(f) 


25(f) 


nn (j-g) 






< 1.5 






■-'■^'-'-0.06 


o q+0.7 

■'•■^-o.s 










< 120 


< 100 


70+*" 
"J-eo 


5200+2600 
ozuu_2200 




630.9/328 


380.3/324 


432.1/324 


572.5/326 


372.3/323 


356.8/322 


360.0/322 



Notes: Results of phenomenological fits with a variety of continuum models. The continuum in Model 1 is assumed to consist of a simple 
powerlaw having a normalisation Nhard and a DISKBB component. The feature around the 6-7 keV range is modelled with a smeared edge 
(Model la) together with a further Gaussian line (Model lb). Model 2 replaces the smeared edge and the Gaussian line with a single relativistic 
line component. M odel 3 replaces the DISKBBan d powerlaw models with the Comptonization code of Poutanen & S vensson ( 1996) and the 
model PEXRIV of iMagdziarz & Zdziarskil jl995t) respectively. For Model 3, a "slab" geometry was assumed and the reflection option in 
compPSwas deactivated. A lower limit for the electron temperature of lOkeV was imposed. The feature is modelled with a broad Gaussian 
(Model 3b) and a relativistic line (Model 3c). Model 3d self consistently convolves the PEXRIV model with the same parameters as the 
relativistic line since they both originate from the same region. All errors are 90 per cent confidence. 



served between the different models, in § 4. The scenario so 
far presented here however, has an edge energy of 5.4 kcV 
which is much less than the value for neutral iron absorption; 
itself a lower limit since iron is likely to be highly ionised in 
the inner parts of the accretion disk. In fact, it is more likely 
that the edge is trying to compensate for a broad emission fea- 
ture (see § 4) and, for this reason, although the model is de- 
tailed in Table 1, we do not believe it to convey any physical 
information. 

Given the obvious presence of a broad line feature in Fig- 
ure 2, we replaced the SMEDGE model with a LAOR line 
profile ( lLaoilll99l1) as expected if emission is coming from 
the inner disk around a black hole. The line energy is con- 
strained to lie between 6.4 and 6.97 kcV, thus encompass- 
ing the full range of possible ionisation states of iron. We 
start in Model 2 with a simple powerlaw emissivity profile 
such that e(r) oc r~^"* . The outer disk radius was frozen 
at the maximum value in the model of 400 rg. This model 
is detailed in Table 1 and resulted in a satisfactory fit with 
(X^/'^ — 432.1/324). Allowing for a broken powerlaw emis- 
sivity profile with indices qin within a radius Vbreak and qout 
beyond, further improved the fit ( A^^ = 13.5 for 2 degrees of 
freedom). This fit has an emissivity index of qin ^ 7 within 
a radius of ^ Qrg and then breaks to qout ^ 2.6. In both 
these instances, the inner radius obtained is very low imply- 
ing that not only is the central object in MAXI J1836-194 a 
stellar mass black hole, it is also likely to be rapidly spinning. 
In all models considered so far, the DISKBB component con- 
sistently required a disk with a temperature of approximately 
0.45 keV. 

It has been suggested by a number of authors that this sim- 
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Fig. 3. — Data/Model ratio for the various phenomenological models de- 
scribed in Table 1. From top to bottom: phabs * smedge * {diskbb + 
powerlaw) with the smeared edge energy frozen at 7.11 keV; phabs * 
{laor + diskbb + powerlaw); phabs * {compps + pexriv); phabs * 
(compps + pexriv + gaussian); phabs * (compps + pexriv + laor); 
phabs * [compps + kdblur (X) pexriv + laor). 

pie view of a multicolour disk is not app ropriate and should 
be rep laced with a broader disk model. iKolehmainen et al] 
( 1201 ll) further suggested that the broad iron line constantly 
seen in a number of stellar mass black hole binaries might be 
an artificial effect caused by the usage of narrower disk com- 
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Inner radius (r ) Emissivily index (q.^) Inner radius (r^) 



Fig. 4. — Contour plots investigating the effect of the inchnation (left) and emissivity index (right) on the inner radius for Model 3d. The centre panel shows the 
contour plot for the inclination versus emissivity index. The 68, 90 and 95 per cent confidence range for two parameters of interest are shown in black, red and 
green respectively. The cross marks the global minima. It is clear that even for the most curved continuum where the emission feature would appear the narrowest 
(i.e. Model 3), the inner extent of the accretion disk, as obtained solely from the phenomenological fit to the iron line profile, is still consistent with not being 
truncated fai' beyond the innermost stable circular orbit for a Schwarzschild black hole (6rg) in the hard/intermediate spectral state observed here. However, due 
to the local minima in these phenomenological models we cannot rule out a disk that is mildly truncated at this stage. 



ponents similar to DISKBB, used here. Despite the noticeable 
problems with this interpretation — iron lines are mostly seen 
in black hole binaries in the low/hard state where the disk tem- 
perature is relatively cold and its contribution to the spectrum 
at ?a 6 kcV is unimportant — we nonetheless test the effect 
this may have on the residuals seen between 4-7 kcV (Fig- 
ure 2), by replacing both the DISKBB, LAOR, and POWER - 
LAW compo nents with compPS (Pou t anen & Svenssonll 19961) 
and PEXRIV (iMagdziarz & Zdziarskil 19951) . The former fully 
characterises the process of Comptonization for a variety 
of coronal geometries, electron distributions and seed pho- 
tons injection geometries. We used the code assuming a 
purely thermal electron distribution in the corona (Gmin=- 
1). The PEXRIV model represents a exponentially cut-off 
power law spectrurrQ reflected from ionized material. When 
the reflection fraction, i? > 0, this models give the sum 
of the illuminating powerlaw hitting the ionised disk to- 
gether with the corresponding reflection component (Fe-edge 
and Compton hump) but it does not include the Fe Ka 
emission line. We initially have the ionisation parame- 
ter = AnF/n crgcms"^, where F is the incident flux 
and n is the number density of hydrogen nuclei), reflection 
fraction and inclination of the PEXRIV component frozen 
at 1000 erg cm s^^, 1 and 25 degrees respe ctively, in line 
with th e values used for XTE J 1752-223 by iNakahira et alj 
(1201 lah . The geometry of compPS is chosen to be a "slab" 
(geom=l) and the covering fraction frozen at 1. Reflection 
from the compPS component is turned off (rel_refl=0). 

This combination resulted in a unsatisfactory fit with 
X^/v = 629.4/327. Allowing the ionization of PEXRIV to 
vary did not resolve the problem (y^ /v = 572.5/326), with a 
broad feature clearly present in the residuals between 4-7 kcV 
(Model 3a in Table 1; Figure 2). Further allowing the reflec- 
tion fraction to change again did not provide a satisfactory 
fit (x^/iy = 546.8/325). The reflection fraction artificially 
shoots up to ^ 2.5 as the PEXRIV model tries to compensate 
for the lack of an Fe-emission line by increasing the depth of 
the iron-edge (see Figure 9). Indeed, adding a broad GAUS- 
SIAN line (Model 3b) resulted in a dramatic improvement, 
with Ax^ = 200.2 for 3 degrees of freedom compared to 
Model 3a. However, a line energy of Ecauss = 5.5 ± 0.1 kcV 
is not consistent with emission of iron. Furthermore, a width 

' The cut-off energy is frozen at 300 keV, similai' to the value used in the 
reflection model REFLIONX described in detail in the next section. 



of ^ 800 cV is highly suggestive of emission from close to 
the black hole where gravitational broadening effects are im- 
portant, again taking us back to the need of a relativistic iron 
line. 

We therefore proceed by replacing the GAUSSIAN line in 
Model 3b with the relativistic line expected around a spin- 
ning black hole. The inclination in the LAOR model is 
frozen at 25 degrees, as with PEXRIV. This model immedi- 
ately improved the quality of the fit (Ax^ = 15.5 for 1 de- 
gree of freedom) and, more importantly, brought the emis- 
sion line energy to a range consistent with emission from iron 
(ELaor ~ 6.6 keV). The inner radius as obtained from the 
iron line profile in Model 3c is again consistent with the cen- 
tral object in MAXI J1836-194 being a rotating black hole. 
Allowing the reflection fraction in the PEXRIV model to be 
free yields R ~ 0.9 ± 0.2, and does not change any of the 
other fit parameters nor does it affect the quality of the fit 
(X^/i^ = 356.4/321). Importantly, Model 3c as it stand is not 
physically consistent. If the Fe line is being emitted in the in- 
ner parts of the accretion disk, as appears to be the case, then 
the other reflection features should also experience the effect 
of strong gravity. For this reason we convolve the PEXRIV 
model — which models the illuminating continuum as well as 
the absorption edge of iron and the Co mpton reflec tion hump 
— with the relativistic kernel KDBLUR jLaod 19911) . We force 
the parameters of KDBLUR to be the same as that of the rela- 
tivistic iron line. This model (Model 3d) result in a slightly 
worst quality of fit Ax'^ = 3.2 for the same number of de- 
grees of freedom as the previous model, however it is now 
physically consistent. Again, allowing the reflection fraction 
in the PEXRIV model to be free resulted in an improvement 
of Ax'^ = 4.4 for 1 degree of freedom, giving R < 0.91 
but did not change the line profile in any way, with the inner 
radius and emissivity index remaining at rin= S.'St'^'^rg and 
Qin = 3.4 ± 0.2, respectively. We also investigated the ro- 
bustness of the fit with respect with the cross-normalization 
constant between the XIS and PIN data which we have so far 
fixed at 1.16. Allowing this to go free barely improved the 
quality of the fit (Ax^ = 1.4 for 1 degree of freedom) and 
recovered a value of 1.14lQ gg, fully consistent with the ex- 
pected value. 

To expedite the computational time, in all incarnations of 
model 3 we have assumed that the inner disk inclination has a 
value of 25 degrees. As a last step in our exploration of these 
phenomenological models, we allow the inclination to be free 
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TABLE 2 

Fits with physically motivated reflection models 



Nn ( xl022cm-2) 

r 

fcTdi,fe(keV) 

■^hard 

Flllum/ FbB 

riff (xlQiS) 

^refbhb 

Nreflionx (XlQ-^) 

Qin 

Qout 

'^break ('"g) 

9 (degrees) 

nn (rg) 
Spin (a) 

X> 



Model 4a 
0.11 ± 0.02 
1.92 ± 0.04 
448+0-006 



< 0.01 

+0.5 
0.4 



4.5 



3.7 ±0.4 

> 5.8 
3.4 ±0.1 



3.6 



+2.3 



2150 
2.54 



0.2 

< 15 

170 
250 
+0.06 
0.20 



379.7/322 



Model 4b 
0.11(f) 
1.92(f) 
0.448(f) 
< 0.01 
4.5(f) 



4 5+0-5 
> 4.4 
3.7 ±0.1 



Model 5 



3.2_ 

< 30 



1+2.1 
■ 0.1 



O+0.8 
^-0.3 



198.5/211 



0.191 



a 

+0.003 



-0.030 
9+0.06 
-0.01 
7+0.02 
-0.01 
041+0-0" 



1.82: 

0.27^ 



1.2 ±0.2 
2.5 
0.31 



+0.7 
1.0 
+0.03 
0.04 



> 6.5 

o 19+0.06 
■^■^^-0.12 
o 4+0.2 

< 23 



9+0.2 
-0.4 



2.2: 

344.9/322 



Model 5b 
0.19 ± 0.03 



1.81 



+0.05 
-0.03 
"+0.02 



0-26lo.03 
030+0-077 



1.2 
2.9 

0.30 



+0.2 
0.3 
1.7 
0.6 
+0.02 
0.04 



3.19 
3.6 



> 7.3 

+0.07 
0.05 
+0.2 

0.1 
< 17 



0.88 ± 0.03 
344.5/322 



Notes; Model 4 is described in XSPEC as phabs * (diskbb + powerlaw + kdblur * reflionx). Model 4b is identical to 4a but we only fit 
the 3-10 keV energy range. Model 5 replaces REFLIONX and DISKBB with the fully self-consistent reflection model REFBHB. In all models 
described so far the kernel from the LAOR line profile was used to account for the gravitational effects close to the black hole. In Model 5b, 
we finally replaces the KDBLUR kernel with the relativistic code RELCONV were the spin is a parameter of the model. In all cases the hard 
emission illuminating the disk is assumed to be a powerlaw with index F. All errors are 90 per cent confidence for one p arameter. " Inner 
radius is not a model parameter and was derived by using the relationship between the spin and ISCO jBardeen et alJI 19720 . It is shown here 
merely to allow for easy comparison. 



and investigate any possible degeneracy this might have on the 
inner accretion radius and emissivity index. Figure |4] shows 
the 68, 90 and 95 percent confidence range for both param- 
eters as a function of inclination as well as the inclination as 
a function of emissivity index. From Fig. 4, it can be seen 
that despite the fact that the inclination is not very well con- 
strained, ranging from anywhere between 5 and 40 degrees at 
the 90 per cent confidence level, the inner radius of the ac- 
cretion disk, as obtained solely from the breadth of the line 
profile using purely phenomenological models, is still consis- 
tent with being at or close to the radius of the Innermost Stable 
Circular Orbit around a rotating black hole. The global min- 
imum which is marked with a black cross in all panels, still 
requires ri„ < 5.6rg at the 90 percent level however we do 
see the presence of a further solution having a radius which is 
marginally consistent with the ISCO of a non-rotating black 
hole together with a very high emissivity index. In Fabian et 
al. (2012), we showed that a single power-law emissivity pro- 
file has only limited validity and that a slope of ^ 3 is a fair 
approximation for an inner disk starting at rjn'^ 2rg. How- 
ever, it severely underestimates the profile within 2rg 
and is therefore a poor probe of the innermost region around a 
rapidly spinning black hole. In that paper, we argued that if q 
= 3 is used, then it will likely yield an upper limit to the inner 
radius and thus a lower limit on the spin if the source is in- 
deed spinning rapidly. The evidence so far points toward the 
black hole in MAXI J1836-194 not being rapidly spinning, 
and therefore, from theoretical arguments we should expect 
an emissivity index close to 3, similar to the value found in the 
the global minima having an inner radius within ri„ < 5.6rg. 
In the following section, we will investigate these possible de- 
generacies fully using a number of self-consistent reflection 
models. 

In this section, we have established beyond any reason- 
able doubt the presence of a strong and mildly broad emis- 



sion line associated with iron in the Suzaku spectra of 
MAXI J1836-194. We have shown that different continuum 
models, as well as phenomenological features (such as the 
blurred edge component) does not eliminate the need for a 
broad emission line. Our efforts to "remove" the need for a 
broad emission line are somewhat artificial since reflection is 
a natural consequence of a system where hard X-rays are im- 
pinging on a cold accretion disk. However, it is only by doing 
so that we can safely rule out any continuum effect on the 
breadth of the iron emission line. In the following section, 
we endeavour to interpret the spectra in a fully consistent and 
physical manner. 

PHYSICALLY SELF-CONSISTENT MODELLING 

In all our previous fits, a broad feature has been shown 
to be robustly present above a thermal-disk and powerlaw- 
like continuum, peaking around the energies expected for 
neutral Fc I emission (~ 6.4 keV) and highly ionised 
H-like FeXXVILya (- 6.97 keV; See Fig. pj. The nat- 
ural explanation for this feature is that it is indeed as- 
sociated with the reprocessing of hard X-ray emission by 
an accretion disk and, as such, the broad feature observed 
is the signature of iron fluorescence that has been rel- 
ativistic broadened due to the strong gravity around the 
central black hole. Similar reflection features are ob- 
served in a wide range of objects ranging from neu- 
tron s tars ("Bhattacharvva & Strohmavei"2007'; 'Cackett et alj 
120081 |2009a b; di Salvo et al. 2009; Reis et al. 2009b), stellar 
mass black holes (Miller 2007; Blumet al. 2009; iReis et all 
i2009a; Hiemstra et al. 2011; Walton et al. 2012) and^AGNs 
fenaka et af.l ll995l; iFabian et al. 2009; Miniutti et aL 2009; 
ISchm oU etai] 120091: TWalton et al..,2010t; iNardini et al.i i20l! ; 



Brenneman et al.ll201 iHReis et al.ll201 lah . 



Due to the immense diagnostic potential of reflection fea- 
tures, a large theoretical effort has been devoted to fully char- 
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Fig. 5. — From top to bottom: Response unfolded vFu spectrum to the 
best-fit Model 5b. The total and blurred REFBHB components are shown as 
cyan and blue solid lines respectively. The powerlaw continuum falls below 
the y-scale; Data/model ratio for: (4a) phabs * (diskbb + powerlaw + 
kdblur ® reflionx); (5a) phabs * {powerlaw + kdblur (g) refbhb) and 
(5b) phabs * {powerlaw + relconv ® refbhb). 



acterisin g the reflection spectrum expected to arise from such 
systems (iLightman & Whit^ 119881: iGeorge & Fabianl 1 1991 ' 



Matt et all 119911: iRo ss & Fa bianI [19931; iZycki et al.l 119941: 
Navakshin et al.ll2000 ; Ballan tvne et al.ll2001l; lRoss & Fabianl 



2005L I2007t iGarcfa & KaUmanI l2010l:TGarcfa et al.1 12011 ? 



Amongst these, the most widely used refl ection model is the 
REFLIONX code of iRoss & Fabianl (120051) . This model self- 
consistently calculates the reflection arising from all ener- 
getically important ionization states and transitions expected 
in disks around black holes. At low ionization parame- 
ters it reproduces the refle ction continuum first described by 
ILightman & Whit3 (|i988!), as well as self-consistently calcu- 
lating the fluorescent lines; at higher ^, lower Z element be- 
comes ionized which results in a softening of the reflection 
spectrum. We start by using a combination of REFLIONX, 
relativistically convolved with KDBLUR, together with disk 
emission (DISKBB) and a powerlaw illuminating continuum. 
The photon indices of the REFLIONX and powerlaw compo- 
nents are assumed to be the same. The outer disk radius is 
assumed to be at 400rg (the maximum allowed by the model) 
and the iron abundance of MAXl J1836-194 is fixed at solar. 
This model give a poor fit with x^/i' ~ 400.6/324. Allowing 
instead for a broken powerlaw emissivity profile results in a 
significant improvement with x^/^ = 379.7/322. This mod- 
els is described in detail in Table 2 (Model 4a) and shown in 
Fig.H 

It is worth noting some similarities and differences between 
the values obtained from the current reflection model and 
those from the phenomenological fits described in the previ- 
ous section. To start, the high ionization parameter found here 
(^ « 2200 erg cm s^^) is similar to that found in the most 
physically motivated version of Model 3 (i.e. Model 3d) and 
suggest an intermediate to high ionization state. The moderate 
to low disk inclination shown in Fig.|4|is also confirmed here, 
where the current reflection fit suggests 9 < 15°. The disk 
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Fig. 6. — Response unfolded uFu spectrum from Model 4b fit in the 3- 
lOkcV energy range. The total, blurred reflection and disk components are 
shown as cyan, blue and red solid lines respectively. The powerlaw contin- 
uum falls below the y-scale. Bottom: Data model ratio to Model 4b. The 
reflection parameters obtained in this range are identical to that obtained by 
modelling the full spectra (see Table 2). 



parameters found here are also in close agreement with those 
found in Model 2, however we cannot compare the parame- 
ters with Model 3 as the disk was modelled assuming a Comp- 
tonization model. A major difference between the fit with the 
reflection over that of the various phenomenological models 
is in the necessity for a broken emissivity profile over that of 
a single powerlaw. In model 4a, the disk extends to within 
^ 2.5rg. The steep emissivity occurs in a very narrow annuli 
between this radius and ^ 4rg, at which point it goes back 
to the value expected from a purely Newtonian geometry. It 
is interesting to note that in the phenomenological model 3d, 
the inner radius obtained with an emissivity index of ^ 3.3 is 
similar to the break radius found here. In the case of the phe- 
nomenological models, the inner radius was obtained from a 
pure emission line . The approach of using a model such as 
the LAOR line profile to obtain spin has many limitations as 
the information imprinted by the effects of strong gravity is 
not limited to an emission line. In fact, such a limitation is 
clear in the inconsistency between the high ionization found 
here (^ ^ 2200 erg cm s^^) as well as that found in Model 3d 
(<^ ^ 5000 erg cm s^^) with the line centroid of ~ 6.48 keV 
found in that same model which indicates lowly ionised iron. 
In Model 3d, a further constraint on the inner radius was in 
place by the act of convolving the PEXRIV model with the 
same kernel as the LAOR line profile. Here, there was an in- 
terplay between broadening of the emission line component 
as well as the iron absorption edge (modelled separately and 
therefore not forced to be physical). The final product was 
the apparent presence of an emission line at ~ 6.48 keV, as 
expected from cold matter with ^ < 100 erg em s^^, together 
with the absorption edge from hot hydrogenic iron Ka FeXXV 
and/or FeXXVl with edge energies > 8.85 keV as expected 
from 500 ^ ^ ^ 5000 erg em s^^. This combination acted to 
model a broad feature with only a mild degree of relativistic 
broadening. REFLIONX on the other hand, by virtue of the 
imposed self-consistency in its atomic physics, modelled the 
same broad feature by allowing a greater degree of broaden- 
ing to be attributed to gravity, and hence resulted in a smaller 
inner radius. 

A question that generally arises when one is dealing with 
reflection models is whether the constraints on the various pa- 
rameters of interest (i.e. inner radius/spin, disk inclination and 
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Fig. 7. — Contour plots investigating the effect of the inclination (top-left), inner (top-centre) and outer (top-right) emissivity indices, PIN/XIS cross- 
normalisation constant (bottom-left), ratio of illuminating to black body flux (bottom-centre) and disk-surface hydrogen density (bottom-right) on the inner 
radius for Model 5a. The latter two being a proxy to the common disk ionisation parameter The 68, 90 and 95 per cent confidence range for two parameters of 
interest are shown in black, red and green respectively. The cross marks the global minima. It is clear that the inner extent of the accretion disk is robust to these 
parameters and is well constrained to ~ 1.9 — 2.3 rg at the 90 per cent level of confidence. 



emissivity profile, etc) are driven by anything other than the 
reflection features. For example; could it be that extreme blur- 
ring, which would suggest a maximally rotating black hole, 
is artificially caused by the model trying to "smooth" the soft 
part of the reflection spectrum to mimic a disk component? To 
test this we ignored the XIS data below 3 kcV where the disk 
emission dominates (see Fig. |2]i and removed the PIN data 
altogether What remained is essentially the broad feature 
which we are associating with iron fluorescence emission. We 
refit Model 4a after freezing the disk temperature and its nor- 
malization as well the neutral hydrogen column density and 
powerlaw index, as these cannot be constrained from the line 
profile alone. We refer to this as Model 4b in Table 2, and 
show this fit in Fig. |6] It is clear that all parameters obtained 
from this narrow energy range are consistent, within errors, to 
that obtained using the full spectra. Allowing F to vary does 
not change this conclusion. 

Until this point, we have made use of REFLIONX which 
is designed to reproduce reflection spectra from the accre- 
tion disks around AGN. However, the disks around stellar 
mass black holes are significantly hotter, resulting in subtle 
differences in the radiation processes, including reprocessing. 
The higher disk temperatures around stellar mass black holes 
means that Compton broadening is of greater importance and 
shoul d therefore be included at the correct level. With this in 
mind, iRoss & FabianI (l2007h developed a modified version of 
the previous grid, REFBHB, in which the atmosphere of the ac- 
cretion disk is is illuminated not only by the hard, powerlaw- 
like corona, but also by a further blackbody radiation intrin- 
sic to the disk. This model self-consistently accounts for the 
disk, relativistic line and reflection continuum present in the 
phenomenological models described in the previous section. 

The parameters of the model are the number density of hy- 
drogen in the illuminated atmosphere, uh, the temperature 



of the accretion disk, kTdisk the index of the (assumed pow- 
erlaw) continuum and the ratio of the total flux illuminat- 
ing the disc to the total blackbody flux emitted by the disc, 
FiUum/ Fbb- The disc reflection spectra is again convolved 
with KDBLUR to account for relativistic effects. This model 
(Model 5a; Table 2 and shown in Figure 5) provides the best 
fit yet to the data with /v = 344.9/322 despite being more 
constrained — compared to purely phenomenological mod- 
els — by virtue of being physically self-consistent. The inner 
accretion disk radius of = 2.2tQ4rg is consistent with 
the values found from the REFLIONX model together with 
a separate disc component. Assuming that this radius is the 
same as the radius of the innerm ost stable circular orb it, we 
can constrain the spin parameter dBardeen et al.|[T972b to be 
a = 0.915t°°^. 

We again investigate the dependence of the inner radius rm 
on a number of key parameters including the inclination an- 
gle, the inner and outer emissivity indices, the cross normal- 
isation between the XIS and PIN data, the ratio of the illu- 
minating powerlaw to the blackbody flux and the hydrogen 
number density at the disk surface. Figure 7 shows that the 
inner radius found here of of Tin = 2.2'^^\rg is extremely 
robust to changes in all the aforementioned parameters. How- 
ever, in order to make a formal constraint on the spin, we re- 
place KDBLUR with the sophis ticated variable-spi n relativis- 
tic smearing model RELCONV (iDauser et al.ll2010l) . Figure 8 
shows that the spin of the black hole in MAXI J1836-194 is 
well constrained to be a = 0.88 ± 0.03 at the 90 per cent 
confidence range. In the following, section we will discuss, 
amongst other things, the current strength and limitation of 
the various models used throughout this work and highlight 
some of the factors contributing to the tight constraint on the 
spin parameter of MAXI J1836-194. 
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Fig. 8. — Goodness-of-fit versus spin for Model 5b. It is clear that the black 
hole in MAXIJ1836-194 is rapidly rotating. A maximally rotating Kerr or a 
static Schwarzschild black hole are both rejected at greater than the 3o" level 
of confidence. 

4. DISCUSSION 

The fractional contribution of the disk component to the to- 
tal 2-20 kcV flux 25 %) places the cu rrent ob s ervati on in 
the hard/intermediate state as defined by iBellonil ( 1201 Oh and 
in an intermediat e state between the low/ hard and high/soft 
state definition of lMcClintock et aTl ( 120061) and. However, we 
can see from figures 5 and 6 that the observed continuum is 
not dominated by the powerlaw-like component expected to 
originate from the corona, but rather it is mostly reflection 
dominated. Similar "reflection-dominated" spectra are seen 
in a number of narrow line S eyfert Is and quasars during their 
low state (e . g. 1H 041 9-577, iFabian et al.l (12005b: NGC 4051, 
iPpnti et al. ("2006^; PG1543-H489. IVignali et alT(E 008): Mrk 
335, Grupeet al, (2008); PG1535-t -547, iBaUaS al. (2008); 
PG21124-059, IScharteletal.1 (l2010h '). In this scenario, the 
majority of the X-ray reprocessing (reflection) occurs in the 
inner region of the accretion disk where strong gravitational 
hght be nding is expected to occur (iMartocchia & Matlill996l 
iMrniutt i & Fabian 2004). Such behaviour is expected as a re- 
sult of strong light bending, where the reflected flux is en- 
hanced over the inner regions as a result of gravitational fo- 
cusing of the X-ray continuum down towards the black hole 
and onto the disk. The decrease in the number of X-rays that 
can escape as part of the continuum thus causes the source to 
appear reflection dominated. 

Wilkins & Fabian (2011) showed that a possible conse- 
quence of strong gravitation effects is an increase in the emis- 
sivity profile of the disk. Classically, the emissivity profile is 
expected to be flat in the region directly below the source, 
while tending to r^^ when r >> h, where the flux re- 
ceived by the disc from the source falls off as the inverse 
square of the distance with a further factor of 1/r arising from 
the cosine of the angle projecting the ray normal to the disc 
plane. However, as brie fly mentioned above and detailed in 
IWilkins & Fabianl(l201 lb . strong gravity can potentially act to 
focus more of the direct continuum into the inner parts of the 
disk as well as increase the disk area being radiated - the latter 
as a consequence of gravitational warping. In this particular 
treatment, these factors causes for a substantial steepening of 
the emissivity profile in the inner regions. 

This is indeed a possible explanation for what is observed in 
MAXI J1836-194, where within a radius of < 4rg the emis- 
sivity index is consistently > 6 (see Table 2) and beyond it 
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Fig. 9. — Difference between the blun'ed REFLIONX component in 
Model 4a (black) to the PEXRIV components in Models 3d (blue) and 
3c (red). The unbluiTed PEXRIV (red) is as expected from a "cold" 
(^ = 70 erg cm s~^) accretion disk being illuminated by a F = 2.18 pow- 
erlaw without being relativistic convolved. Compare this to the component 
expected from a moderately ionised (g = 5200 erg cm s^^; blue) disk at 
3.3rg illuminated by a F = 2.09 powerlaw and the REFLIONX model 
having both the absorption edge as well as the iron emission line (black). 

goes closer to the classical value of 3. Both the reflection- 
dominated spectrum and the high emissivity profile seen here 
suggest that the primary X-ray continuum is located within 
a few gravitational radii of the black hole. Looking at Fig- 
ure 1, it is indeed possible that the corona briefly "collapsed" 
down close to the black hole causing the decrease in the hard 
X-ray flux seen during this Suzaku observation. If this is the 
case, state transitions (at least that between the high/soft and 
low/hard) should be seen to be much more tightly associated 
with changes in the corona as opposed to physical changes in 
the accretion disk. 

In this work we have investigated a variety of possible mod- 
els striving towards a physically motivated interpretation for 
the observed spectrum. It is worth stressing that equally good 
fits - and at times (albeit not in this work) statistically better 
fits - can be obtained with a purely phenomenological com- 
bination of components severely lacking in physical consis- 
tency. A case in point is Model 3c where the presence of 
the - clearly broad - iron line requires the emission to oc- 
cur from deep within 6rg, yet the direct X-ray and reflection 
continuum appear somehow exempt from the effect of strong 
gravity. When we account for relativistic effects, the ioniza- 
tion parameter of the disk increases by nearly two orders of 
magnitude from ^ w 70 to « 5000 erg cm s~^. This change 
is accompanied by a hardening of the photon index F possi- 
bly due to the powerlaw trying to compensate for the stronger 
soft emission and the weaker Compton hump of the blurred 
ionised reflector in comparison to the unblurred, cold reflec- 
tor (see Fig.|9l). 

A further point to note is that the apparent decrease in the 
equivalent width of the relativistic line profile (from 270 cVto 
180 cV from model 3c to model 3d) is clearly a consequence 
of the, physically inconsistent, lack of gravitational blurring. 
From Fig.|9]we can see that the absorption edge in PEXRIV 
becomes much more smooth and symmetric after the compo- 
nent is convolved with KDBLUR and the ionization increases. 
Due to the decoupling between the LAOR line profile and 
PEXRIV, this smooth edge can act somewhat like an emission 
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line conspiring against the LAOR line component and thus de- 
creasing the equivalent width of latter In order to couple the 
emission line strength with the absorption edge and all other 
reflection features, we replaced the phenomenological com- 
bination of LAOR + PEXRIV with the REFLIONX reflection 
grid. The first thing to notice from Table 2 (model 4) is that 
this recovered the high ionization value from Model 3d but 
now in a self-consistent manner. Nonetheless, the strong re- 
quirement for an emission line with an equivalent width of 
at least ?» 180 cV (Model 3d) is highly indicative of a broad 
line, adhering to the strong criterion of Reis et al. (20fOI) ar- 
guing against an accretion disk truncated far from the ISCO. 

The precise value of the neutral hydrogen column density 
towards MAXl J1836-194 is not known, however, based on 
Swift/XRT data which extends to much lowe r energies as 
compa red to the current Suzaku observation, iKennea et al.l 
(1201 Ih showed that a likely value is (2 ± 0.4) x 10^^ cm"^. 
A similar value is found here using our best model. Figure 9 
highlights how differing models impact the energies below 
^ 1 kcV where we would begin to see the curvatures expected 
from such a low column density. The fact that our data cuts 
off at 1.2 kcV means that, in order to break the degeneracies 
between the - essentially disk - continuum (see Fig. 2) and 
neutral column density, one must use a self consistent model 
for the full spectra. Such is the case for REFBHB. 

It is well known that in the hot inner regions of an ac- 
cretion disk the observed disk spectrum suffers from the 
effects of electron scattering which results in an observed 
(colour) temperature, Tcoi, which is higher than the effec- 
tive blackbody temp erature, Thh, by approximately a factor 
of = Tcoi/Tbb dRoss et al.lll992l) . This "colour correc- 
tion factor" has been show n to have a value of 1.7 ± 0.2 
dSh imura &Takaharall995h for a wide range in luminosity, as 
long as the disk effective temperature remains below ^ 1 keV 
(IPavis et al. 2005), as is the case here. Above this tempera- 
ture, disk self ionisation can lead to an increas e in fr,ni how- 
ever i t is found to be consistently below ^ 3 (iMerloni et al.l 
120001) . The effective temperature in the REFBHB (Model 5) 
of fcTs = 0.27j^Q Q^ kcV is precisely as expected from the 
value of the colour temperature from the DISKBB component 
in Model 4 {kT^/kT^ = 1.7 ± 0.1) . 

Lastly, although we cannot directly compare ^ between 
REFLIONX and REFBHB, we can estimate this value based 
on the various parameters output in Table 2. Since the flux 
of the blackbody in the model is related to its temperature 
by the F^b oc relation, we find Fbb = (S.S^'^J g x 
10^^ erg cm"2 s"^ Combining this F/num /Fbb = 1.2±0.2 
gives an illuminating flux of (6.5l^'4) x 10^^ erg cm~^ s^^. 
The ionization parameter, which is defined as ^ = 
AnF/n crgcms~^, is then found by using n = nn = 
2.5t? o X lO^^H cm-3. In fliis manner, we find ^ = 3300 ± 
1500 erg cm s^^, in perfect agreement with the values found 
for both REFLIONX in both Models 4a and 4b and the bluiTed 
PEXRIV in Model 3d. 

All current methods of measuring black hole spin rely on 
both the assumption that the accretion disk extends to the in- 



nermost stable circular orbit, as well as that emission within 
this radius is negligible. The latter is indeed valid for the 
standard model of black hole accretion (IShakura & SunyaevI 
Il973h . where within the ISCO - in the region often referred 
to as the plunging region - there is no angular momentum 
transport, and the region cannot support an X-ray corona to 
irradiate the material that is ballistica lly plunging o nto the 
black hole. Ho wever, it w as shown by iKrolid (119991) and in- 
dependently by iGamml^ (17999) that when magnetic fields are 
considered, the B-fields within the ISCO may be amplified to 
a point where the magnetic energy density can be comparable 
to the rest mass energy of the accreting material and may lead 
to the creation of an active inner X-ray corona. IShafee et al.l 
( 120081) investigated the effect of magnetic torque within the 
plunging region around a non-rotating black hole and con- 
cluded that "...magnetic coupling across the ISCO is relatively 
unimportant for geometrically-thin discs". 

A further study specifically aimed at addressing the robust- 
ness of the iron line/reflection fitt ing technique in diagnosin g 
black hole spin was presented by [Reynolds & FabianI (120081) . 
The authors used a high resolution three-dimensional magne- 
tohydrodynamic simulation of a geometrically-thin accretion 
disc to show that the density of the plunging material drops 
precipitously over a very small radius within the ISCO. This 
sudden drop in density results in the material being highly 
photoionised and suppresses any significant iron line emission 
as well as all other reflection features from within the ISCO. 
The study by Reynolds & Fabian (2008) concluded that for 
a non-rotating black hole where the ISCO is at 6rg, the re- 
flection edge - defined by the authors as the innermost ra- 
dius from which significant reflection emission is seen - is 
at approximately 5.8rg. Furthermore, the discrepancy be- 
tween the true ISCO and the inferred radius diminishes as 
one considers more rapidly rotating black holes as is the case 
in MAXI J1836-194 (see Fig. 8). From an observational per- 
spective, it is also worth noting that strong support for the 
presence of an "inner edge" in the accretion disk surrounding 
black ho les is provided by d ecades of empirical evidence as 
shown in lSteiner et al.l ( 120101) . 

The culmination of this work is that the recently discovered 
system, MAXI J1836-194, is indeed a stellar mass black hole 
binary, having a central black hole rotating with a spin pa- 
rameter of a = 0.88 ± 0.03 (90% confidence). This strong 
constraint is a result of being able to successfully, and most 
importantly self-consistently, model the reflection features 
clearly present in the Suzaku spectra. 
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